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This paper discusses the application of the CLS to generate probabilistic 
loads for use in the PSAM nonlinear evaluation of stochastic structures under 
stress (NESSUS) finite element code. The CLS approach allows for quantifying 
loads as mean values and distributions around a central value rather than 
maximum or enveloped values typically used in deterministic analysis. NESSUS 
uses these loads to determine mean and perturbation responses. These results 
are probabilistically evaluated with the distributional information from CLS 
using a fast probabilistic integration (FPI) technique to define response 
distributions. The main example discussed describes a method of obtaining 
probabilistic (dependent) pressure, temperature and centrifugal steady state 
load descriptions and stress response of the second-stage turbine blade of the 
Space Shuttle main engine (SSME) high-pressure fuel turbopump (HPFTP) . 
Additional information is presented on the on-going analysis of the high 
pressure oxidizer turbopump discharge duct (HPOTP) where probabilistic dynamic 
loads have been generated and are in the process of being used for dynamic 
analysis. Example comparisons of load analysis and engine data are furnished 
for partial verification and/or justification for the methodology. 

Figure 1 depicts the component loading analysis.^ The description of the 

dependent random variables are obtained from a multilevel physical model that 
includes an engine system influence model and component load models of the 
turbine to describe the local conditions on the turbine blade. The 
independent loads are quantified at the appropriate model — engine or 
component - to best define the basic variables in the system, as well as 
maintain the required correlation between the dependent loads on the turbine 
blade. Both duty cycle and related probabilistic variations are accounted for 
in these models. The specific turbine blade model. Figure 2, is used both for 
part of the load modeling-thermal and pressure distribution - as well as the 
structural analysis. The structural model combines the loads with a 
probabilistic description of geometry, material property, and material 
orientations. The work herein discussed is an extension of the analysis 
performed under the PSAM contract where probabilistic loads were not 

considered . 

Figure 3 is a cross section of the HPFTP turbine showing the overall 
configuration and the hot gas coolant network. The hot gas flow enters from 
the fuel preburner at 11 and discharges at 15. The second stage blade and 

rotor are between 14 and 15. Two geometric variations play a key role in the 

distribution of the coolant flow around the turbine blade, the interstage seal 
at 10 and the aft platform seal at 8. 

The HPFTP turbine blade loads are dependent to the engine system level major 
component (turbine) dependent loads. The major component dependent loads are 
in turn dependent to a set of engine independent loads, including the engine 
operating parameter loads and the engine hardware parameter loads. A summary 
and categorization of these loads are listed in figure 4. The independent 
system level loads and hardware variations are in the first category listed. 
The major component dependent loads that are calculated in the engine model 
are the next level. The component level dependent loads are pressure, 
temperature and centrifugal. Additional local independent loads from the 
seals are internal to the turbine. 

*Work performed under NASA Contract NAS3-24382 
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Probabilistic methods are built into these models and are described in the CLS 
discussion in this document. For this analysis, the Gaussian moment method 
was employed to illustrate the application of the probabilistic methodology 
The pressure, temperature and centrifugal loads are correlated with higher 
level system independent loads. Theoretically, one could obtain correlation 
fields between pairs of the loads and then decompose them as part of the PSAM 
solution scheme. These fields are difficult to obtain and decompose. The 

marginal distribution method is used as a better solution to the load 

description and interface to the PSAM codes. Figure 5. Figure 6 outlines the 
process of generating the HPFTP turbine component loads marginal distribution. 

At this time, validation of these loads is limited to the standard 

instrumentation available from engine tests. Pump speed and turbine discharge 
temperatures are regularly recorded as well as the independent loads at the 
pump interfaces. Figure 7 has typical plots of test data from a series of 
tests that can be used to estimate 2 sigma variations of independent loads. 
Figure 8 has typical dependent pump speed variations. A sample from a flight 
is also furnished of the 3 pump speeds at 104% power level. The variation in 
speeds are caused by a combination of hardware and duty cycle variations of 
primarily the low pressure oxidizer pump inlet pressure and the resulting 

cavitation of the HPOTP. Table 1 furnishes a correlation of calculated and 
measured variations of turbine speed and discharge temperature for both high 
pressure turbopumps. There is excellent correlation with speed and slightly 
poorer correlation with temperature. Part of the temperature error is 
attributed to instrumentation error that is not readily accounted for. 

Table 2 summarizes the set of random variables, how they affect the structural 
model and their mean and standard deviation values. Figure 9 shows two points 
on the structural model where cumulative distributions of effective stress 
were calculated using NESSUS and FPI codes. 

PSAM is performing a series of verification analysis on SSME hardware. CLS is 
being used to develop the loads for these analyses to furnish realistic input 
to PSAM and to demonstrate the integrated use of CLS and PSAM. Vibration 
loads for the HPOTP discharge duct. Figure 10, have been completed. This duct 
is attached to the HPOTP at two points. Zone G, and at the main injector (MI) 
at Zone A of the engine. The MI is part of a combustor that generates random 
vibrations. The pumps generate both random and sinusoidal loads. These 
vibration loads are dependent on power level of the engine and pump power and 
speed. The structure of the engine couples the responses of the components 
such that sinusoidals from the high pressure pumps are found throughout the 
engine at attenuated levels compared to the point where they are generated. 
Table 3 summarizes the 38 individual load components that are being applied to 
the duct. Figure 11 shows one of the six normalized PSD shapes used for the 
random loading. For this analysis, the random load composite has been modeled 
as a normalized shape that is scaled by G£ ms . The sinusoids are depicted 

In Figure 12. The frequencies generated by a pump are correlated to the 
synchronous speed, IN. The IN, 2N and 3N response levels are typically caused 
by rotor unbalance and rubbing or nonlinear effects. The 4N is caused by 
fluid loads on the four primary blades of the inducer. The variations are 
based on engine test data. The sinusoid amplitude mean levels are correlated 
with pump speed, but the standard deviation magnitudes are uncorrelated. The 
sinusoid level distribution is modeled as lognormal, see Figure 13. The 
variations of the sinusoidal magnitudes have a coefficient of variation of .5 
to over 1.0. The high variations are roughly equally caused by component 
hardware variations and test to test variations. Figure 14 is a recent study 
of random composite levels coefficients of variations for several environments 
on the SSME. This information will be used in the development of generic 
dynamic physical load models. 
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FIGURE 1. COMPONENT LOADS 



FIGURE 2. COMPONENT MODELS FOR LOADS AND STRUCTURAL ANALYSIS 
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FIGURE 3. SSME HP FTP TURBINE HOT GAS AND COOLANT FLOW NETWORK 
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FIGURE 4. COMPOSITE LOADS OF THE HPFTP TURBINE BLADE 
FOR THE PROBABILISTIC STRUCTURAL ANALYSIS 
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PURPOSE: 




• PREPARING COMPOSITE LOADS FOR STRUCTURAL ANALYSIS. 

• MARGINAL DISTRIBUTION LOADS USED DIRECTLY AS INPUT TO PROBABILISTIC STRUCTURAL 
ANALYSIS CODES, E.G. NESSUS & FPI . 

• ADVANTAGES: 

• MUCH SIMPLER THAN CORRELATION FIELD APPROACH. 

• CORRELATIONS TO THE ENGINE INDEPENDENT LOAD LEVEL AS REQUIRED.. 

• DIRECTLY USABLE FOR PERTURBATION LOADING IN PSAM. 

• DEFINITION: 

• MARGINAL DISTRIBUTION FOR A COMPONENT LOAD IS A DISTRIBUTION FOR THE COMPONENT 
LOADS WITH ONLY ONE INDEPENDENT LOAD VARYING AND THE REST OF THE INDEPENDENT 
LOADS STAYING CONSTANT. 

• FOR A LINEAR SYSTEM 
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riGURE 5. MARGINAL DISTRIBUTION METHOD 
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FIGURE 6. MARGINAL LOAD PROCESS 
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FIGURE 7. TYPICAL INDEPENDENT LOAD VARIATIONS 
FOR A PHASE II SSME SERIES OF TESTS 
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TABLE I. CALCULATED VS. MEASURED VARIATIONS IN PUMP SPEED AND TEMPERATURE 
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NO 

RANDOM VARIABLE 
DESCRIPTION 

| TYPE 

1 

1 feh | 

| QUANTITIES | 
! AFFECTED I 

MEAN 

STANDARD 

DEVIATION 



1 

MATERIAL AXIS 

1 


-0.087266 

0.067544 




ABOUT Z 

1 

| MATERIAL | 

RADIANS 

RADIANS 





| MATERIAL 






2 

MATERIAL AXIS 

i 

]ORIENTATION| 

-0.034907 

0.067544 




ABOUT Y 

I AXIS 


RADIANS 

RADIANS 





i 

| ANGLES I 





3 

MATERIAL AXIS 

I VARIATIONS 


+0.052360 

0.067544 




ABOUT X 

i 


RADIANS 

RADIANS 



4 

ELASTIC MODULUS 

i 


T8.38E6 

0.4595E6 





| ELASTIC 


KSI 

KSI 





i 

| ELASTIC i 





5 

POISSON’S RATIO 

| PROPERTY 


0.386 

0.00965 





1 

| CONSTANTS | 





6 

SHEAR MODULUS 

[ VARIATIONS 


18.63E6 

0.46575E6 





i 


KSI 

...... KSI 



7 

GEOMETRIC LEAN 

i 


0.0 

0.14 




ANGLE ABOUT X 

i 



DEGREES 





| GEOMETRY 

| NODAL I 





8 

GEOMETRIC TILT 

i 


0.0 

0.14 




ANGLE ABOUT Y 

I VARIATIONS 

(COORDINATES! 


DEGREES 



9 

GEOMETRIC TWIST 

i 

i 


0.0 

0.30 




ANGLE ABOUT Z 

i 



DEGREES 



TO 

MIXTURE RATIO 

i 







LIQUID HYDROGEN/ 

| INDEPENDENT 


6.00 

0.02 




LIQUID OXYGEN 

1 








| LOAD 






11 

FUEL INLET 

1 


30.00 psia 

5.00 




PRESSURE 

| VARIATIONS 
1 

| PRESSURE | 





12 

OXIDIZER INLET 

1 

| DEPENDENT 

1 - 1 

100.00 psia 

26.00 




PRESSURE 

i 

| TEMPERATURE | 







| LOADS ARE 






13 

FUEL INLET _ 

1 


38.5* R 

0.5 




TEMPERATURE 

| TURBINE BLADE 
1 

| CENTRIFUGAL | 





14 

OXIDIZER INLET 

| PRESSURE 


167.0* R 

1 .33 




TEMPERATURE 

1 

| LOAO | 







| TEMPERATURE 






15 

HPFP EFFICIENCY 

1 


1.00 

0.008 





| AND SPEED 






16 

HPFP HEAD 

1 


1.0237 

0.008 




COEFFICIENT 

1 


, 




17 

COOLANT SEAL 

| LOCAL 


1.00 

0.1 




LEAKAGE FACTOR 









| GEOMETRY 

(TEMPERATURE | 





18 

HOT GAS SEAL 

1 


1.0 

0.5 




LEAKAGE FACTOR 

I FACTORS 


1 




TABLE 2. LIST OF 2nd STAGE BLADE RANDOM VARIABLES 
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FIGURE 9. STRESS RESPONSE FOR HPFTP TURBINE BLADE 
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HPOTP 



FIGURE 10. HPOTP DISCHARGE DUCT 
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SEQUENTIAL 

NO. 

1 

_J_ 

RANDOM VARIABLE 
DESCRIPTION 


I MEAN | 

ST. D. 

1 

| 

COV 

1) 

1 

ZONE G - X AXIS, 

PSD 

| 222.0 | 

163.0 

| 

0.73 

2) 

1 

ZONE G - Y AXIS, 

PSD 

I 73.5 | 

59.4 

| 

0.808 

3) 

1 

ZONE G - Z AXIS, 

PSD 

1 73.5 | 

59.4 

| 

0.808 

4) 

1 

ZONE A - X AXIS, 

PSD 

1 22.5 | 

4.5 

| 

0.20 

5) 

1 

ZONE A - Y AXIS, 

PSD 

| 54.0 | 

10.8 

| 

0.20 

6) 

1 

ZONE A - Z AXIS, 

PSD 

I 69.5 | 

13.9 

| 

0.2 

7) 

1 

OXIDIZER PUMP SPEED 

I 2940.53 | 

41 .167 

| 

0.014 

8) 

1 

FUEL PUMP SPEED 


| 3707.08 | 

37.07 

| 

0.01 

9) 

1 

1 

1 

DAMPING 

ZONE A OXIDIZER TURBOPUMP 
OXIDIZER PUMP SINE AMPLITUDES 

| 0.033 | 

1 1 

1 | 

.005 

1 

1 

| 

0.15 

10) 

1 

X DIRECTION 

IN 

1 .30 | 

.120 

| 

0.4 

ID 

1 


2N 

I 0.30 | 

.045 

| 

0.15 

12) 

1 


4N 

1 1.5 | 

.45 

I 

0.3 

13) 

1 

Y DIRECTION 

IN 

I 0.60 | 

0.30 

| 

0.5 

14) 

1 


2N 

I 0.70 | 

0.28 

| 

.40 

15) 

1 


4N 

1 2.6 | 

.78 

| 

.3 

16) 

1 

Z DIRECTION 

IN 

1 0.5 | 

.225 

I 

0.45 

17) 

1 


2N 

I 0.70 | 

.140 

1 

0.20 

18) 

1 

1 

4N 

FUEL PUMP SINE AMPLITUDES 

I 0.70 | 

1 | 

.140 

1 

| 

.20 

19) 

1 

X DIRECTION 

IN 

1 .35 | 

1.050 

| 

.3 

20) 

1 

Y DIRECTION 

IN 

I 0.80 | 

0.280 

| 

0.35 

21) 

1 

1 

1 

1 

Z DIRECTION IN 

ZONE G - MAIN INJECTOR 
OXIDIZER TURBOPUMP 
SINE AMPLITUDES 

1 1.20 | 
1 1 

1 1 

1 | 

0.96 

1 

1 

1 

| 

0.3 

22) 

1 

X DIRECTION 

IN 

I 1.35 | 

1 .35 

| 

1 .0 

23) 

1 


2N 

I 1.50 | 

0.75 

I 

0.5 

24) 

1 


3N 

1 1.10 | 

0.495 

I 

0.45 

25) 

1 


4N 

1 n.o | 

2.75 

| 

0.25 

26, 27) 

1 

X-Y DIRECTION 

IN 

1 1.9 | 

1 .71 

| 

0.9 

28, 29) 

1 


2N 

1 1.6 | 

0.96 

I 

0.6 

30, 31) 

1 


3N 

I 0.75 | 

.225 

| 

0.3 

32, 33) 

1 

1 

4N 

FUEL PUMP SINE AMPLITUDES 

1 5.5 | 

1 | 

3.30 

1 

| 

0.6 

34) 

1 

X DIRECTION 

IN 

I 0.65 | 

0.2275 

| 

0.35 

35, 36) 

1 

X-Y DIRECTION 

IN 

I 0.45 | 

1 .35 

| 

.3 

37. 38) 

_L_ 

Y-Z DIRECTION 

2N 

1 0.45 1 

0.180 

J 

.4 


Note: 1) P.S.D. units are in G 2 

2) Pump Speed units are in Radians/SCE 

3) Sinusoidal amplitude units are in GRMS 


TABLE 3. BASIC RANDOM VARIABLES AND THEIR STATISTICS (INPUT) 
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FIGURE 12. RANDOM AND SINUSOIDAL VIBRATION 
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FIGURE 13. HPOTP SYNCHRONOUS MAGNITUDE 
DISTRIBUTION AT CONSTANT POWER 



FIGURE 14. SSME ZONAL VIBRATION COMPOSITE GRMS 
COEFFICIENTS OF VARIATION 







